Nonlinear Two-Dimensional Green's Function in Smectics 



E A. Brcner 1 , V. I. Marchcnko 2,3 and D. Pilipenko 1 , 
1 Institut fur Festkorpcrforschung, Forschungszcntrum Jiilich, D-52425 Jiilich, Germany 

2 Kapitza Institute for Physical Problems, Russian Academy of Sciences, ul. Kosygina 2, Moscow, 119334 Russia 

3 Moscow Institute of Physics and Technology, Institutskii per. 9, Dolgoprudnyi, Moscow region, 141700 Russia 

(Dated: March 23, 2010) 

The problem of the strain of smectics subjected to a force distributed over a line in the basal 
plane has been solved. 



The asymptotic expressions for strains around isolated 
defects in smectics at long distances are characterized by 
the exponent a |. If a < 0, the linear theory is ap- 
plicable. If a = 0, (edge dislocation Greens func- 
tion [H), the linear solution is valid at small action am- 
plitudes, whereas nonlinear effects become important at 
larger amplitudes. If a > 0, nonlinear effects should be 
taken into account even for extremely weak actions. In 
this paper, we report the solution of the nonlinear prob- 
lem of the two-dimensional Greens function {a = 1/2) 
[H]. Let us consider a smectic sample with the thick- 
ness L that is sandwiched between solid undcformable 
walls parallel to smectic layers (see figure). A force uni- 
formly distributed along the y axis (normal to the figure 
plane) with the linear density F is applied at the cen- 
ter (x = z = 0) of the smectic sample. The energy of 
the small strains of the smectic sample is given by the 
expression (see Eq. (44.13) in (H[) 
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where u is the displacement of the layers along the smec- 
tic axis z, A is the elastic modulus, and A is the micro- 
scopic length parameter. In our problem, the maximum 
displacement uq is reached in the force application line. 
In terms of the new function / = u/uq and new coordi- 
nates z = z/L and x = x/L, where e = uq/L, Eq. (fT]) is 
represented in the form 
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where (3 = (X/eL) 2 , a = f - (f) 2 /2, and the dot and 
prime mean differentiations with respect to z and x re- 
spectively. Thus, it is necessary to determine the func- 
tion / that provides the minimum of energy (J5|) , is equal 
to unity at x = z = 0, and is equal to zero at the edges 
of the smectic layer z = ±1/2. The force is given by the 
expression F = F+ — F- , where 



a zz dS = ^ALe 3/2 / adx 
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In the macroscopic problem, curvature oc /3 can be 
neglected even for the case of a negligibly small force 




FIG. 1: Strain of the smectic layer subjected to the lin- 
early distributed force(e = 0.1)applied at the thick point. The 
dashed line is the boundary of the region, where a noticeable 
strain appears. The thin line is the boundary of the compres- 
sion region(cr < 0). 



F ~ XA^J X/L, where the amplitude uq is larger than 
the distance between smectic layers ~A<L. It is in- 
teresting that integrals ((2]) and ((3]) can be calculated in 
this case even without the complete solution, because 
the strain field is divided into two regions. In the first 
region, where the material is compressed, the problem 
can be solved analytically. In the second region, owing 
to the Helfrish instability (see the problem in Section 
44 in [H[), rotary states almost without stresses appear 
instead of tension. Here, noticeable stresses exist only 
inside microscopically thin twin boundaries Q. 

The figure shows the strain pattern in the smectic layer 
obtained by numerically solving the problem. In the com- 
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pression region, the function / is a quadratic function of 
the variable x: f(x,z) = fo(z) + f2(z)x 2 . It is easy to 
verify that the equilibrium equation a — (/'&)' = cor- 
responding to the extremum of energy @ has such an 
exact solution satisfying the necessary conditions. It can 
be represented in the parametric form 
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The compression region a < corresponds to the range 
— x TO < x < x m , at z > 0, where y/7r. The 

thin line in the hgure is the boundary taking into account 
the displacement of smectic layers, where a is zero. Un- 
der this line, down to the dashed line, strains result in 
the appearance of two twins @ in which inhomogeneous 
strains increase in the direction to the force application 
level. Here, stresses are small and are likely caused by 
the finiteness of the grid used for the calculation. 

When the tension of the twin boundaries is disregarded 
( at P = 0) @, F_ = 0. Then, F = F+ = (8/37r)ALe 3 / 2 . 
Correspondingly, 



u 



2/3 



l}' 3 , a 



F 
X 



At z <C 1, according to |QJ, u = u + Su, where 
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, 5u <C u . 
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Let us assume that the displacement near the force ap- 
plication point has the form Su = — z 1 / 3 a 2 / 3 ?/'(u), where 
v = xaT 1 / 3, z~ 2 / 3 . Note that Suis independent of L. In 
this case, the equilibrium equation is the ordinary differ- 
ential equation 



((2V> - ivtp' + 3(V>') 2 ) (3^ - 2v))' = 0, 
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where the prime means differentiation with respect to v, 
and has a trivial first integral. The integration constant 
is zero, because the expression in the second parentheses 
in Eq. ([7]) is zero for solution ((6]). The expression in the 
first parentheses is proportional to stress a. It is conve- 
nient to represent the solution of the equation a = in 
the parametric form if) = (4 + t 3 )/2t, and v = (1 + t 3 )/t 2 . 
The intervals -1 < t < and < t < 2 1 / 3 -> z > cor- 
respond to the regions z < and z > respectively. 
The integration constant is chosen from the condition 



tp (3/2 2 / 3 ) = 3/2 1 / 3 at t = 2 1 / 3 (matching with solution 
© on the line x m (z) ~ 3(az 2 /4) 1 / 3 ). Near the z = line 
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At the x = line at z < (|u| ->■ 0) 
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If a <C A, onlinear asymptotic expressions are valid at 
|z| 3> A 3 /a 2 , and |x| ^> A 2 /a, and the Greens function 
of the linear approximation is applicable at smaller dis- 
tances [l|. If a 3> A, small-strain approximation ([T|) is 
violated at \z\ ~ |x| ~ a. 

Stresses in the problem under consideration exist only 
in the smectic compression region and rotary adjust- 
ment in a certain bounded region occurs instead of ten- 
sion. Such a character of nonlinear response implies that 
the compression field and, correspondingly, the uo value 
(with the change L — > 2L+), as well as the asymptotic 
expression for Su, remain unchanged in the general case, 
where the force is applied not to the center, but at any 
distance L + in the direction of the force action from the 
undeformable wall. The boundary conditions on the op- 
posite side and, generally speaking, beyond the compres- 
sion region affect only the adjustment structure at dis- 
tances of about the sample sizes. This work was sup- 
ported by the German-Israeli Foundation and the Rus- 
sian Foundation for Basic Research (project no. 09-02- 
00483). 
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